基于 尺度 不 变 特征 点 匹配 的 太阳 磁场 图 像 配 准 与 定位 方法 ' 


杨 盼 '， 曾 电光"”*， 刘 锁 ?， 郑 胜 1， 林 钢 华 ?， 何 慧 灵 1 
(4. 三 峡 大 学 理学 院 ， 湖 北 宜昌 443000 ; 2. 中 国 科学 院 太 阳 活 动 重点 实验 室 ， 
北京 100012) 


摘要 不 同 太 阳 观 测 站 间 存 在 观测 时 间 、 观测 范围 、 观测 设备 等 差异 ， 为 了 方便 科学 研 
究 ， 有 必要 找到 不 同 观测 站 间 太 阳 图 像 自动 配 准 与 定位 的 方法 。 提 出 了 基于 尺度 不 变 特征 点 
匹配 的 太阳 磁场 图 像 配 准 与 定位 方法 。 首先 ， 对 原始 图 像 进行 对 比 度 增强 、 降 采样 等 预 处 理 ; 
其 次 ， 采 用 尺度 不 变 特征 检测 算法 提取 两 图 的 尺度 不 变 特征 点 ， 然 后 ， 寻 找 同名 点 的 方法 
- 对 两 幅 图 像 特征 点 进行 粗 配 准 和 粗 定 位 ， 最 后 ， 基 于 粗 定 位 区 域 实现 大 阳 磁场 图 像 的 精确 
- 配 准 和 定位 ,开展 了 不 同时 段 磁 场 图 像 的 配 准 和 定位 实验 ， 对 匹配 点 对 数 、 匹 配 准确 率 、 匹 配 
误差 等 参数 进行 了 定量 分 析 ， 实 验 结果 表明 ， 该 方法 可 以 自动 、 准确 、 快速 地 实现 太阳 磁场 
图 像 的 配 准 与 定位 。 
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6 发 活动 ， 须 监测 从 太阳 光 球 到 日 园 的 磁场 和 流 场 结构 及 其 演化 过 程 "。 目 前 ， 国 内 外 已 经 奸 
一 有 多 个 太阳 磁场 观测 设备 。 如 国内 怀柔 太阳 观测 基地 (Huairou Solar Observing 
c Station,HSOS) 的 口径 35cm 太阳 磁场 望远镜 (Solar Magnetic Field Telescope, SMFT) P, 
国外 太阳 动力 学 天 文 台 C Solar Dynamics Observatory, SDO ) 的 日 震 及 磁场 成 像 仪 
(Helioseismic and Magnetic Imager, HMI) B 等 。 不 同 太阳 观测 站 间 存 在 观测 时 间 、 观 测 范围 、 
观测 设备 等 差异 .为 了 方便 科学 研究 ， 提 高 太阳 物理 数据 使 用 效率 ， 有 必要 找到 不 同 观测 站 
间 太 阳 图 像 自动 配 准 的 方法 。 另 外 ， 当 两 观测 设备 的 观测 范围 不 同时 ， 还 需要 实现 “小 图 
像 ”在 “大 图 像 ” 中 的 自动 定位 例如， 怀柔 太阳 观测 基地 的 太阳 磁场 望远镜 得 到 的 是 局 部 
太阳 磁场 图 像 ， 然 而 ， 太 阳 动 力学 天 文 台 的 日 震 及 磁场 成 像 仪 得 到 的 是 全 日 面 太阳 磁场 图 
像 。 前 者 的 观测 范围 远 远 小 于 后 者 ， 这 就 需要 找到 一 种 能 够 将 局 部 太阳 图 像 在 全 日 面 图 像 中 
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快速 定位 的 方法 。 
目前 已 有 一 些 关 于 太阳 图 像 配 准 方面 的 研究 工作 ， 如 文 [4] 提 出 了 基于 修正 矩 的 亚 像素 
AGA EV RUE EXE, SCISPRAASE T fri BS SIFT 算法 对 太阳 图 像 进行 了 配 准 。 这 些 都 是 针 
对 同一 观测 站 的 序列 图 像 进行 的 ， 其 成 像 设备 位 置 固定 ， 图 像 内 容 在 短 时 间 内 变化 较 小 ， 
位 移 、 旋 转 和 仿 射 变换 较 少 。 不 同 太阳 观测 站 的 太阳 图 像 由 于 在 观测 时 间 、 观 测 范围 、 观 测 设 
等 方面 存在 较 大 差异 ， 需 要 适应 力 更 强 的 图 像 配 准 方 法 进行 配 准 。 已 有 的 图 像 配 准 方法 很 
多 ， 如 基于 灰 度 的 配 准 方法 Wl 基于 特征 的 配 准 方法 忆 等 等 。 基 于 灰 度 的 配 准 方法 具有 精度 
高 的 优点 ， 但 计算 复杂 度 高 ， 且 对 目标 的 旋转 、 形 变 、 遮 挡 以 及 灰 度 变化 等 很 敏感 。 基 于 尺度 
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适应 能 力 。 太 阳 磁 场 图 像 含有 许多 内 部 结构 ， 特 征 点 较 明 显 ,本 文 提出 了 基于 尺度 不 变 特征 


点 匹配 的 太阳 磁场 图 像 配 准 与 定位 方法 ， 开 展 了 怀柔 太阳 观测 基地 的 局 部 太阳 磁场 图 像 在 
太阳 动力 学 天 文 台 日 震 及 磁场 成 像 仪 的 全 日 面 太阳 磁场 图 像 中 自动 配 准 与 定位 的 研究 ， 并 
做 了 检测 实验 证 明 方 法 的 有 效 性 。 


1. 算法 原理 


= 尺度 不 变 特征 点 匹配 算法 对 局 部 太阳 图 像 在 全 日 面 太阳 中 自动 配 准 与 定位 ， 主 要 包括 
= 尺度 不 变 特征 提取 算法 的 特征 点 检测 、 特 征 点 初步 匹配 、 特 征 点 精确 匹配 、 图 像 变换 矩阵 求解 
= 4 个 步骤 。 
四 1.1 SIFT 特征 点 检测 

1.1.1 尺度 空间 极 值 检测 
- 尺度 空间 理论 是 检测 不 变 特征 的 基础 。 文 [12] 证 明了 高 斯 卷 积 核 是 实现 尺度 变换 的 唯一 
= 变换 核 。 
E 一 幅 二 维 图 像 在 不 同 尺度 下 的 尺度 空间 可 以 表示 为 图 像 与 高 斯 核 卷 积 : 


L(x, y,0) = G(x, y,o)*I(x, y) ， (1) 


其 中 ，(x;y) 为 图 像 点 的 像素 坐标 ，Z(x,y) 为 图 像 灰 度 值 ，G(X,y,0) 为 高 斯 核 函 数 。 
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G(x, y,0o)= e 


200° 


(2) 
其 中 ，c 为 尺度 空间 因子 ， 它 是 高 斯 正 态 分 布 的 方差， 反映 了 图 像 被 平滑 的 程度 ; 
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L(x, y,0) 为 图 像 的 尺度 空间 。 


为 高 效 地 在 尺度 空间 检测 稳定 的 极 值 ， 文 [10] 使 用 尺度 空间 中 差分 高 斯 DoG 的 极 值 作 
为 判断 依据 。DoG 算 子 定义 如 下 : 


D(x,y,0)=[G(x,y,ko)— G(x,y,0)]* I(x, y) 
=L(x,y,ko)— L(x, y,o) 


(3) 
其 中 , 为 相 邻 两 个 尺度 空间 的 比例 因子 。 

DoG 局 部 极 值 的 检测 是 将 每 个 像素 的 DoG 值 跟 同一 尺度 的 周围 邻 域 8 个 像素 和 相 邻 尺 
度 对 应 位 置 的 周围 邻 域 9x2 个 像素 总 共 26 个 像素 进行 比较 。 仅 当 被 检测 点 的 DoG EKF 26 
个 像素 点 或 小 于 26 个 像素 点 时 ， 才 将 其 判断 为 极 值 点 并 保存 ， 以 进行 后 续 计算 。 

1.1.2 确定 关键 点 位 置 及 尺度 

通过 拟 合 三 维 二 次 函数 以 精确 确定 关键 点 的 位 置 和 尺度 ， 从 而 可 得 到 原 图 像 尺 度 不 变 
特征 提取 算法 的 候选 点 集合 .其 中 ， 对 比 度 低 的 点 对 噪声 比较 敏感 ， 而 位 于 边缘 上 的 点 难以 
准确 定位 。 为 了 确保 尺度 不 变 特征 提取 算法 特征 点 的 稳定 性 ， 须 从 尺度 不 变 特征 提取 算法 候 
选 点 集合 将 这 两 类 点 剔除 。 

1.1.3 关键 点 方向 确定 

通过 确定 关键 点 的 方向 ， 可 以 使 特征 描述 符 以 与 方向 相关 的 方式 构造 ， 从 而 使 算 子 具 
旋转 不 变性 。 关键 点 的 方向 利用 邻 域 像素 的 梯度 分 布 特性 确定 。 对 每 个 高 斯 图 像 ， 每 个 点 


L(xy) 的 梯度 方向 0(x,y) 可 以 由 下 式 计算 9; 


E L(x, y +1)- L(x, y -1) 
A(x, y)= ta 人 1, y) | 7 


1.1.4 特征 向 量 生成 
为 了 确保 描述 符 的 旋转 不 变性 ， 首 先 需 旋转 坐标 ， 保 持 坐 标 方向 与 极 值 点 方向 一 致 ， 
然后 将 平面 内 可 能 分 布 的 360" 梯度 方向 划 为 8 个 方向 范围 ， 每 个 方向 范围 各 包含 45° 。 在 子 
区 域 为 4x4 的 情况 下 ， 共 计 有 4x4x8=128 个 数据 ， 即 生成 128 维 描 述 符 。 方 向 与 信息 联合 不 
仅 增强 了 抗 品 性 ， 还 消除 了 尺度 变化 、 旋 转 和 变形 的 影响 * 呈 。 
1.2 特征 点 初步 匹配 
完成 特征 点 检测 后 ， 可 利用 寻找 同名 点 的 方法 对 两 幅 图 像 特征 点 进行 粗 匹 配 5 I, PABX 
式 距 离 4 作 为 图 像 间 的 相似 性 度量 以 确定 同名 点 。 欧 式 距 离 二 维 表示 为 


d= TE L2! i (5) 
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首先 在 待 配 准 图 像 中 找 出 两 个 与 参考 图 像 中 对 应 某 极 值 点 的 欧 氏 距离 最 短 的 点 ; 然后 
以 两 个 值 中 小 的 值 作为 分 子 ， 大 的 值 作为 分 母 ， 若 该 分 数值 在 给 定 的 闵 值 范围 内 ， 那 么 建 
立 对 应 点 的 同名 点 关系 ， 称 为 匹配 点 。 阔 值 的 选取 将 直接 影响 同名 点 的 数量 。 立 值 越 大 ， 获 
得 的 匹配 点 越 多 ， 可 能 的 误 匹 配 也 越 多 ; 六 值 过 小 时 ， 匹 配点 会 很 少 ,。 大 量 实验 表明 ， 立 值 
在 0.4—0.6 Zz [H] fg FED, 
1.3 特征 点 精确 匹配 及 变换 矩阵 求解 

在 完成 特征 点 粗 匹配 后 ， 可 利用 随机 抽样 一 致 性 算法 巴 消 除 错 误 的 匹配 点 ， 从 而 提高 
配 准 精 度 。 随 机 抽样 一 致 性 算法 ， 可 以 在 一 组 包含 “外 点 ”的 数据 集中 ， 采 用 不 断 欠 代 的 方 
法 ， 寻 找 最 优 参数 模型 ， 不 符合 最 优 模 型 的 点 ， 被 定义 为 “外 点 ”。 该 算法 最 早 由 文 [9] 提 
出 马 ， 在 图 像 配 准 以 及 拼接 上 得 到 广泛 应 用 。 

一 般 地 ， 两 幅 图 像 间 的 变换 可 以 由 变换 矩阵 表示 吕 1: 


x h x 
v= h, h, h; y 
1 h, h 1 
(6) 


其 中 ， Cx', yO M y) 分 别 为 两 幅 图 像 中 对 应 点 的 像素 坐标 ，h〔 二 0,1,2...7) 为 变换 


和 矩阵 的 8 个 未 知 系数 。 

随机 抽样 一 致 性 算法 可 对 特征 点 进行 精确 匹配 并 求 变换 矩阵 。 首先 随 机 选择 4 组 对 应 点 
组 成 一 个 随机 样本 并 计算 变换 抢 阵 五 ,对 假设 的 每 组 对 应 计算 距离 g， 计 算 与 瓦 一 致 的 内 
点 数 ， 选 择 有 最 多 内 点 数 的 石 ， 在 数目 相等 时 ， 选 择 内 点 标准 方差 最 小 的 那个 解 中 ,最 后 根 
据 图 像 间 变换 矩阵 石 ， 可 以 计算 出 局 部 太阳 图 像 的 4 个 顶点 及 中 心 点 在 全 日 面 图 像 中 的 像 
素 坐 标 ， 从 而 完成 局 部 太阳 图 像 在 全 日 面 太 阳 图 像 中 的 配 准 与 定位 。 


2. 算法 设计 及 实验 分 析 


2.1 算法 设计 
基于 尺度 不 变 特征 匹配 的 太阳 磁场 图 像 配 准 与 定位 方法 流程 如 下 : 
C1) 对 图 像 进行 预 处 理 ， 包 括 降 采 样 、 黑 白 反 相 、 对 比 度 增强 等 ， 对 图 像 降 采样 到 
1024*1024 左右 ; 
(2) 通过 尺度 不 变 特征 提取 算法 提取 图 像 特征 点 ; 
G) 采用 寻找 同名 点 的 方法 对 两 幅 图 像 特征 点 进行 粗 匹 配 ; 
CA) 利用 随机 抽样 一 致 性 算法 去 除 错误 的 匹配 点 ， 并 估计 两 幅 图 像 之 间 的 最 佳 变 换算 
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(5) 根据 仿 射 变换 和 矩阵， 实现 局 部 太阳 图 像 在 全 日 面 太 阳 图 像 中 的 初步 定位 ; 
C6) 利用 初步 定位 的 图 像 区 域 和 怀柔 局 部 图 像 进行 精确 配 准 。 
算法 的 具体 流程 如 图 1。 


1 太阳 图 像 配 准 与 定位 流程 图 


Fig .1 Flow chart of solar image registration and positioning 


2.2 实验 结果 与 分 析 
2.2.1 不 同 观 测 站 太阳 磁场 图 像 的 配 准 与 定位 

以 怀柔 太阳 观测 基地 的 局 部 太阳 磁场 图 像 在 SDO/HMI 的 全 日 面 太阳 磁场 图 像 中 的 自 
动 匹配 与 定位 展开 实验 。 图 2 给 出 了 原始 的 全 日 面 太阳 图 像 及 局 部 太阳 图 像 。 其 中 ， 图 2(a) 
给 出 日 震 及 磁场 成 像 仪 外 得 到 的 全 日 面 太阳 磁场 图 像 ， 拍 摄 时 间 2012 年 11 月 12 日 05 时 48 
分 00 秒 。 图 2(b) 给 出 35cm 太阳 磁场 望远镜 外 得 到 的 局 部 太阳 磁场 图 像 ， 拍 摄 时 间 2012 
年 11 月 12 日 05 时 55 分 34 秒 。 两 者 的 拍摄 时 间 非 常 接 近 。 值 得 指出 的 是 ， 两 者 的 观测 范围 
和 分 辨 率 有 很 大 差异 。 怀 柔 太阳 观测 基地 的 局 部 太阳 图 像 仅仅 是 全 日 面 太阳 图 像 中 很 小 一 的 
部 分 。 全 日 面 图 像 的 尺度 是 4096x4096， 怀 柔 站 局 部 图 像 的 尺度 830x992, 


(a) 全 日 


ELA ER REG LER s (b) 局 部 太阳 磁场 图 像 
图 2 原始 的 太阳 磁场 图 像 


Fig. 2 Original solar magnetic field image(Full-disk Solar magnetic field image in(a)and local solar magnetic 


field image (in(b))) 


图 2 中 两 幅 图 像 的 对 比 度 不 够 。 在 进行 后 续 的 图 像 配 准 前 ， 有 必要 通过 对 比 度 增 强 等 对 


201805.00414v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


两 幅 图 像 进行 预 处 理 。 另 外 ， 由 于 日 震 及 磁场 成 像 仪 的 全 日 面 图 像 分 辩 率 太 高 ， 为 了 提高 配 
准 速度 ， 先 做 降 采 样 处 理 ， 将 尺度 从 4096x4096 降低 到 1024x1024。 最 后 ， 部 分 日 震 及 磁场 
成 像 仪 的 全 日 面 太 阳 磁 场 图 像 的 正 负极 性 与 怀柔 观测 站 的 相反 。 因 此 要 做 黑白 反 相处 理 。 经 
过 降 采 样 、 黑 白 反 相 及 对 比 度 增强 处 理 后 的 全 日 面 太阳 磁场 图 像 如 图 3 Ca) 。 经 过 对 比 度 增 
强 后 的 局 部 太阳 图 像 如 图 3 (b) 。 对 比 图 2 (b) 和 图 3 CbO 可 以 发 现 ， 预 处 理 后 的 太阳 图 
像 比 原始 图 像 更 加 清晰 ， 这 有 利于 后 续 的 特征 提取 与 图 像 配 准 。 


[b] 


(@) 全 日 面 太阳 磁场 图 像 (0) 局 部 太阳 磁场 图 像 
图 3 经 过 预 处 理 后 的 太阳 磁场 图 像 


Fig.3 Full-disk solar magnetic field image(in(a))and local solar magnetic field image (in(b))after pre-treatment 


随后 ， 采 用 尺度 不 变 特征 提取 算法 提取 图 像 特征 点 。 图 4 给 出 了 经 过 预 处 理 后 的 局 部 太 
阳 图 像 的 尺度 不 变 特征 检测 结果 。 一 共 检 测 到 了 204 个 尺度 不 变 特征 点 。 由 图 可 见 ， 这 些 特 
征 点 较 好 地 代表 了 局 部 太阳 图 像 的 特征 。 
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图 


4 经 过 预 处 理 后 的 


局 部 太阳 图 


像 的 尺度 不 变 特征 检测 


Fig.4 Scale-invariant feature detection of localized solar image after pre-treatment 


在 对 图 3 中 的 两 幅 图 像 进行 尺度 不 变 特 和 
的 全 日 面 图像 一 共 检 测 到 了 3447 个 尺度 不 变 特征 点 ， 局 部 太阳 图 像 一 共 检 测 到 了 204 AR 
度 不 变 特征 点 。 接 下 来 ,采用 寻找 同名 点 的 方法 以 及 随机 抽样 一 致 性 算法 找 出 两 幅 图 中 相互 


配 准 的 特征 点 对 。 如 图 5， 


FE 检 测 后 ， 需 要 将 两 幅 图 的 特征 点 配 准 。 图 3 中 


找到 了 20 个 相互 配 准 的 特征 点 对 《图 中 用 红色 直线 连接 ) 。 


根据 这 些 配 准 特征 点 对 ， 可 初步 估计 出 两 幅 图 像 之 间 的 变换 矩阵 A, HAER: 


$9.1375 0.0583 727.0 
H, = fb .0573 -0.1303 675.9 . (7) 


T S dE BE Hn] A H 


0 


0 


局 部 图 像 的 4 个 顶点 及 中 心 点 在 全 日 面 图 像 中 的 像素 位 置 


1 


~ 


从 而 完成 了 局 部 图 像 在 全 日 面 图 像 中 的 粗 定位 (图 5)。 


图 5 局 部 图 像 在 


AHE 


i 图 像 中 的 配 准 与 定位 


Fig.5 Registration and positioning of local image in full-disk solar image 
对 日 震 及 磁场 成 像 仪 全 
降低 匹配 与 定位 精度 。 为 此 ， 将 日 震 及 磁场 成 像 仪 全 日 面 原始 图 像 中 绿色 四 边 形 对 应 的 最 小 


外 接 和 矩形 区 域 ( 见 图 5) 与 怀柔 局 部 太阳 图 像 做 精确 匹配 ， 即 可 确定 局 部 太阳 图 像 在 全 日 


面 图 像 中 的 精确 范围 。 


(图 中 用 红色 直线 连接 ) 。 根 据 这 些 配 ; 准 特 行 


中， 具体 如 (8) 式 ; 


精确 


日 面 图 像 进行 降 采 样 处 理 可 以 大 大 提高 计算 速度 ， 但 同时 也 会 


匹配 后 的 结果 如 图 6， 一 共 找 到 了 25 个 相互 配 准 的 特征 点 对 


点 对， 可 估计 两 幅 图 像 之 间 的 最 佳 变换 矩阵 
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20.5425. 0.2270 547.2349 
H, = Boon -0.5300 662.1596 . (8) 
0 0 1 


由 变换 矩阵 肪 可 以 看 出 ， 两 图 像 间 主要 存在 缩放 、 旋转 、 平 移 等 几何 变换 ， 不 存在 透射 
变换 ， 且 可 以 计算 局 部 图 像 的 4 个 顶点 在 全 日 面 图 像 中 的 精确 位 置 分 别 为 (2909,2704) ~ 
(2371,2484) 、 (2560,2045) ~ (3097,2265) (图 6 中 用 绿色 四 边 形 标 出 ) ， 从 而 完成 了 
局 部 图 像 在 全 日 面 图 像 中 的 精确 定位 另外， 还 可 以 计算 局 部 太阳 图 像 相对 于 全 日 面 图 像 的 
旋转 角度 约 为 157.9?( 顺 时 针 〉; 水 平方 向 放大 比例 为 1.71 、 垂 直方 向 放大 比例 为 1.74。 


6 局 部 图 像 在 全 日 面 图 像 中 的 精确 配 准 与 定位 


Fig.6 Accurate registration and positioning of local image in full-disk solar image 


7 给 出 了 局 部 放大 后 的 全 日 面 太阳 图 像 (图 7 Ca) ) 与 经 过 变换 和 矩阵 用 变换 后 的 局 
部 太阳 图 像 (图 7 b) ) 。 由 图 可 见 ， 两 图 的 形状 、 角 度 等 基本 一 致 ， 但 细节 方面 略 有 差异 。 
两 图 形状 、 角 度 等 基本 一 致 证 明了 本 文 算 法 的 准确 性 。 两 图 细节 方面 略 有 差异 是 因为 两 者 是 
由 不 同 观测 站 得 到 的 ， 在 观测 时 间 、 望 远 镜 参 数 等 方面 存在 差异 所 致 。 但 这 些 细微 的 差异 不 
会 影响 局 部 太阳 图 像 在 全 日 面 图像 中 的 配 准 与 定位 。 
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(a) 局 部 放大 后 的 全 日 面 太 阳 磁 场 图 像 ; (b) 变 换 后 的 局 部 太阳 磁场 图 像 
图 7 全 日 面 太阳 磁场 图 像 与 经 过 变换 后 的 局 部 太阳 磁场 图 像 
Fig.7 Full-disk solar magnetic field image ( in( a ) )and transformed local solar magnetic field image( in ( b ) ) 


2.2.2 不 同时 间 段 的 匹配 结果 

为 了 验证 算法 的 有 效 性 ， 对 71 组 不 同时 段 的 太阳 图 像 进行 配 准 ， 并 给 出 相应 的 配 准 点 
对 数 与 正确 配 准 点 对 数 ， 结 果 如 图 8。 其 中 ， 红 色 实 线 圆圈 为 实际 的 配 准 点 对 数 ， 绿 色 虚 线 
点 为 正确 配 准 的 点 对 数 ， 蓝 色 虚 线 五 角 星 为 配 准 的 正确 率 ( 即 正确 配 准 的 点 对 数 与 实际 配 准 
的 点 对 数 之 比 )， 全 日 面 图 像 分 辨 率 均 为 1024x1024。 由 图 8 可 以 看 出 ， 大 部 分 情况 下 正确 配 
准 的 点 对 数 与 实际 的 配 准 点 对 数 相 等 ， 即 配 准 准确 率 为 100%; 少数 情况 下 正确 配 准 的 点 对 
数 小 于 实际 的 配 准 点 对 数 ， 即 存在 一 定 的 错误 配 准 。 出 现 错误 配 准 的 原因 在 于 : 怀柔 局 部 太 
阳 图 像 和 日 震 及 磁场 成 像 仪 全 日 面 太阳 图 像 ， 一 个 是 地 面 拍摄 的 ， 一 个 是 空间 拍摄 的 ,在 大 
气 沸 流 等 干扰 下 ， 怀 柔 局 部 太阳 图 像 会 变 得 模糊 ，SIFT 算法 能 检测 到 的 特征 点 减少 ， 进 而 
导致 配 准 点 对 数 和 配 准 准 确 率 的 下 降 。 图 9 给 出 了 局 部 太阳 图 像 较 模糊 时 的 配 准 情况 ， 当 怀 
柔 局 部 太阳 图 像 较为 模糊 时 ， 配 准点 对 数 与 配 准 准确 率 较 低 。 
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8 不 同时 间 太阳 磁场 图 像 的 配 准 情况 


Fig.8 The registration of the solar magnetic field image at different time 
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局 部 太阳 磁场 图 像 较 模 糊 时 的 配 准 情况 


Fig.9 Registration results of a blurry local solar magnetic field image 


2.2.3 仿 射 变换 对 太阳 磁场 图 像 配 准 结果 的 影响 

不 同 观测 站 的 太阳 磁场 图 像 间 除了 常见 的 缩放 \ 旋转 和 平移 变换 ， 很 多 情况 下 还 存在 比 
较 明 显 的 仿 射 变换 。 本 节 研 究 仿 射 变换 对 太阳 磁场 图 像 配 准 结果 的 影响 ， 通 过 定量 分 析 并 计 
算 配 准 实验 的 误差 。 从 原始 日 震 及 磁场 成 像 仪 全 日 面 图 像 中 “图 2(a)) 抠 出 一 块 矩形 子 图 像 。 
然后 ， 将 矩形 子 图 像 进行 一 定 的 仿 射 变换 ， 最 后 ， 将 仿 射 变换 后 的 子 图 像 来 与 降 采 样 后 的 


SAMAK (图 3(a)) 做 配 准 ， 将 变换 矩阵 计算 得 到 的 子 图 像 在 全 日 面 图 像 中 的 位 置 与 子 
图 像 在 全 日 面 图 像 中 的 已 知 位 置 相 减 ， 即 可 准确 地 计算 配 准 误差 。 配 准 误差 了 定义 为 


4 
YG =F) FHV) 
n= ， (9) 


4 


Ep, Ges Yo) AO o y ) 分 别 为 子 图 像 的 4 个 顶点 在 全 日 面 图 像 中 的 已 知 位 置 和 经 过 变换 


和 矩阵 预测 的 位 置 。 

仿 射 变 换 包 含 平移 、 旋 转 、 缩 放 、 错 切 等 变换 。 图 10 给 出 了 缩放 为 0.6 倍 的 情况 下 ， 配 准 
误差 (单位 : 像素 ) 随 错 切 系数 和 旋转 角度 的 变化 情况 。 由 图 10 可 以 看 出 ， 当 错 切 系数 较 小 时 ， 
配 准 误差 较 小 ， 而 当 错 切 系 数 增 大 到 一 定 程度 如 超过 0.3 后 ， 配 准 误差 将 迅速 增 大 ， 旋 转 
角度 对 配 准 结果 影响 较 小 。 图 11 给 出 了 旋转 角度 为 45 度 的 情况 下 ， 配 准 误差 随 错 切 系数 和 
缩放 系数 的 变化 情况 。 由 图 11 可 以 看 出 ， 缩 小 或 放大 2 倍 以 内 ， 缩 放 对 配 准 误差 的 影响 较 
小 ， 超 出 此 范围 后 配 准 误差 将 明显 增 大 。 这 是 因为 缩小 变换 需要 丢弃 一 些 图 像 信 息 ， 而 放大 
变换 则 需要 插值 来 增 大 图 像 。 当 缩小 导致 丢弃 的 信息 太 多 或 放大 插值 的 信息 太 多 时 会 对 特征 
点 的 检测 与 配 准 造成 影响 。 综 上 所 述 ， 本 文采 用 的 SIFT 配 准 算法 对 旋转 、 缩放 、 错 切 变 换 等 
均 具 有 一 定 的 鲁 棒 性 .为 了 更 好 地 反映 仿 射 变换 对 配 准 精度 的 影响 ， 本 节 的 结果 都 是 通过 初 
步 定位 得 到 的 (对 应 算法 步骤 a~e) ， 通 过 精确 定位 可 以 大 大 减 小 配 准 误差 。 
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图 10 配 准 误差 随 错 切 系数 和 旋转 角度 的 变化 


Fig.10 The change of registration error varies with the shear factor and the rotation angle 
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图 11 配 准 误差 随 错 切 系 数 和 缩放 系数 的 变化 


Fig.11 The change of registration error varies with the shear factor and zoom factor 


3. 结论 


zl 


采用 基于 尺度 不 变 特征 点 匹配 的 方法 对 怀柔 观测 站 的 局 部 太阳 磁场 图 像 在 太阳 动力 学 
天 文 台 日 震 及 磁场 成 像 仪 的 全 日 面 图 像 中 的 配 准 与 定位 问题 进行 了 研究 ， 结 果 表明 ， 该 算 
法 对 光照 \ 旋 转 、 尺 度 缩放 等 有 很 强 的 鲁 棒 性 ， 对 不 同 观测 站 得 到 的 共有 不 同 观测 条 件 \ 不 同 
观测 尺度 等 的 太阳 磁场 图 像 的 配 准 有 效 且 提高 了 太阳 物理 数据 使 用 效率 。 
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Images Based on Scale-invariant Feature Point Matching 
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Abstract : There are differences among different solar observation stations in terms of 
observation time, scope and equipment. For the convenience of scientific research, it is necessary 
to work out a way to automatically match and locate the solar images of different observation 
stations. A method to match and locate solar magnetic field image based on scale-invariant feature 
point matching is put forward. Firstly, pre-treatment for the original image is carried out by 
contrast enhancing and down sampling; Secondly, the feature scale-invariant feature points are 
extracted using scale-invariant feature point detection algorithm; Next, the feature points of the 
two images are matched and located roughly by seeking their homologous points; Finally, ,Based 
on coarse located area ,the solar magnetic image can be matchend and located exactly . The 
registration and positioning experiments of the solar magnetic field image in different time periods 
is carried out, and number of matching point pairs,matching logarithm, matching accuracy and 
matching error are defined to quantitatively analyze the experimental results. As proved by 
experiment results, the presented method can match and locate the solar magnetic images 
automatically, accurately and quickly. 

Keywords:Solar image; Image registration; Automatic positioning; Scale-invariant; Feature 


extraction 


